Optimal N-to-M Cloning of Quantum Coherent States 
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The cloning of continuous quantum variables is analyzed based on the concept of Gaussian cloning 
machines, i.e., transformations that yield copies that are Gaussian mixtures centered on the state 
to be copied. The optimality of Gaussian cloning machines that transform N identical input states 
into M output states is investigated, and bounds on the fidelity of the process are derived via a 
connection with quantum estimation theory. In particular, the optimal N-to-M cloning fidelity for 
coherent states is found to be equal to MN/{MN + M — N). 
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Cloning denotes an operation by which the unknown 
state of a system is copied. For reasons rooted at the lin- 
earity of quantum mechanics, it turns out that when the 
system to be copied is quantum mechanical, cloning can- 
not be performed exactly Then, a natural question 
that arises is "to what extent can the copies resemble 
the original, in accordance with quantum mechanics?" 
pL This is the problem of optimal quantum cloning, 
which has now been extensively studied for quantum bits 
l^-Q , and, more generally, for d-level systems ]8|-|lfj|] . The 
present paper investigates the question of optimal cloning 
for continuous quantum variables. Examples of contin- 
uous variables include the position and momentum of 
a particle, or the two quadratures of a quantized elec- 
tromagnetic field. We will consider a quantum system 
described in terms of two canonically conjugate opera- 
tors with continuous spectra (referred to as x and p, with 
eigenvalues x and p, respectively). Precisely because they 
are conjugate, x and p cannot be both copied exactly. 
Nevertheless, approximate cloning can be achieved if the 
copies are not required to be exact. Then, the issue of op- 
timal cloning amounts to find the best tradeoff between 
position and momentum errors induced by cloning. 

In this paper, we shall consider N — > M symmet- 
ric Gaussian doners (SGC), defined as a linear com- 
pletely positive map Cn,m transforming N identical 
replicas of an unknown quantum state belonging to 
an infinite-dimensional Hilbert space TL into M > N im- 
perfect clones. The joint state of these clones pu = 
CV,A/(|^® iV )(?/'® Ar |) is required to be supported on the 
symmetric subspace of 7i® M , and is such that the par- 
tial trace over all outputs except one is the bi-variate 
Gaussian mixture 
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where the integral is performed over all values of (3 — 
(x + ip) I V2 in the complex plane (K = 1), and the oper- 
ator D((3) = exp(/3at — j3*a) achieves a displacement of x 



in position and p in momentum, with a and a) denoting 
the destruction and creation operators, respectively [ pT| . 
Thus, the copies yielded by a SGC are affected by an 
equal Gaussian noise <r 2 = er 2 = a 2 N M on the conjugate 
variables x and p. (It will turn out that the resulting 
cloning fidelity / = {ij}\p\\ijj) is invariant for all coherent 
states of x and p.) The symmetry of the doner also ob- 
viously implies that the M copies are characterized each 
by the same density operator p\. 

The issue of the duplication (N = 1, M = 2) of quan- 
tum information carried by a continuous variable has 
been treated in a previous paper [^2[ , where an explicit 
Gaussian 1 — > 2 cloning transformation was proposed. It 
was shown that the noise variance induced by this doner 
is <7 2 2 = 1/2, so that the resulting cloning fidelity for 
coherent states is fx t % = 2/3. This fidelity is invariant 
under translations and rotations in phase space, so that 
this Gaussian doner can be thought of as the analogue 
for coherent states of the universal cloning machine for 
quantum bits [||. The present work investigates the op- 
timality of this 1 — ► 2 doner, and extends these consider- 
ations to N — > M continuous doners. More specifically, 
we address the question of "how close" the output state 
[Eq. ([j])] can be from the input state We find that 
a lower bound on the noise variance cr 2 ^ M is given by 
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implying in turn that the optimal N 
for coherent states is bounded by 



MN 



N,M = 



MN +M -N 



(2) 

M cloning fidelity 
(3) 



First, let us demonstrate that the bound (g) is achieved 
with the 1 — > 2 SGC derived in (l2), so that the latter 
is optimal for coherent states (optimality was only con- 
jectured in |12|). Our proof is directly connected to the 
problem of simultaneously measuring a pair of conjugate 
observables on a single quantum system. It is known 
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(see e.g. |13 ) that any attempt to measure x and p si- 
multaneously on a quantum system is constrained by the 
inequality 
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where a 2 (N) and a 2 (N) denote the variance of the mea- 
sured values of x and p, respectively, when N replicas of 
the state are available. (The case where N > 1 will be 
considered later on.) So, the best possible simultaneous 
measurement of x and p with a same precision satisfies 
cr 2 (l) = CTp(l) = 1. Compared with the intrinsic noise 
of a minimum- uncertainty wave packet a 2 — az = 1/2, 
we see that the joint measurement of x and p effects an 
additional noise of minimum variance 1/2 g|. Now, let 
a coherent state \a) be processed by a 1 — » 2 SGC, and 
let x be measured at one output of the doner while p is 
measured at the other output. As cloning should obey 
inequality m), we must have 



Ax 2 Ap 2 > 1 
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where Ax 2 (Ap 2 ) refers to the usual variance of observ- 
able x (p) measured on p\. Using Eq. (TP), it gives 



(5x 2 +a 2 )(6p 2 
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where Sx 2 (Sp 2 ) is the intrinsic variance of x (p) mea- 
sured on the input state, while o~\ 2 is the noise variance 
induced by the doner. Now, using the uncertainty prin- 
ciple 5x 2 5p 2 > 1/4 and the identity a 2 + b 2 > 2Va 2 b 2 , 
we conclude that the noise variance is constrained by 



A,2 > K-2 = 1/2 
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Let us now estimate a 2 N x , that is, the variance of an 
optimal joint measurement of x and p on N replicas of a 
system. From quantum estimation theory p5[ , we know 
that the variance of the measured values of x and p on 
a single system, respectively <r 2 (l) and cr 2 (l), are con- 
strained by 

g x (j 2 x (l) + g P cTp(l) > g x Sx 2 + g p 6p 2 + ^Jg x g p (9) 

for all values of the constants g x ,g p > 0. Note that, 
for each value of g x and g p , a specific POVM based 
on a resolution of identity in terms of squeezed states 
(whose squeezing parameter r is a function of g x and 
g p ) achieves this bound (see [fl5|| ). Moreover, when mea- 
surement is performed on ./V independent and identi- 



is reduced by a factor 
. So, applying N times 



cal systems, the r. h. s. of ( 
TV -1 , as in classical statistics [ 
the optimal single-system POVM is the best joint mea- 
surement when N replicas are available since it yields 
a 2 x {N) = AT-V^l) and a 2 p {N) = N' 1 ^!). Hence, us- 
ing Eq. ([|) for a coherent state (Sx 2 — 5p 2 — 1/2) and 
requiring a 2 (N) = a 2 (N) , the tightest bound is obtained 
for g x = Op. It yields ~o 2 N x = 1/N, which, combined 
with Eq. (@), gives the minimum noise variance induced 
by cloning, Eq. (§). 

It is now easy to compute the fidelity of the optimal 
N — > M SGC when a coherent state |a) is copied. Using 
Eq. (0) and the identity |(a|a')| 2 = exp(— \a — a'\ 2 ), we 
obtain 
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implying that the doner presented in |l^| is optimal. 

Let us now consider the general problem of optimal 
N — > M Gaussian cloning. Our proof is connected to 
quantum state estimation theory similarly to what was 
done for quantum bits in |14j|, the key idea being that 
cloning should not be a way of circumventing the noise 
limitation encountered in any measuring process. More 
specifically, our bound relies on the fact that cascading 
a N — » M doner with a M — > L doner results in a 
N — > L doner which cannot be better that the optimal 
N — > L doner. We make use of the property that cas- 
cading two SGCs results in a single SGC whose variance 
is simply the sum of the variances of the two component 
SGCs (see Appendix) . Hence, the variance ~a 2 N L of the 
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In particular, if the M 
L — > oo, 



L SGC must satisfy a NL < a N M 



L doner is itself optimal and 
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Since the limit of Cn.m with M — ► oo corresponds to 
a measurement ||, Eq. @ implies that cloning the N 
replicas of a system before measuring the M resulting 
clones does not provide a mean to enhance the accuracy 
of a direct measurement of the N replicas. 



which results in Eq. (||) . As expected, all coherent states 
are copied with a same fidelity. (Note, however, that 
this property does not extend to all states of H.) Equa- 
tions (||) and (||) are consistent with the known result 
for a 1 — > 2 continuous doner, i. e., a 2 2 = 1/2 and 
/i,2 = 2/3 |12|. In addition, they yield the obvious result 



"n n = an d /n,n = 1, confirming that the optimal 
N — > N cloning map is just the identity. Furthermore, 
they fulfill the natural requirement that the cloning fi- 
delity increases with the number of input replicas. For 
instance, considering a kN — » kM SGC with a positive 
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integer k, we find that — < (and 



a-k m>°)- At 

the limit N — > oo, we have / — > 1, VM, that is, clas- 
sical copying is allowed. Finally, for M — > oo, that is, 
for an optimal measurement, we get / — > N/(N + 1). 
In particular, it implies that the best simultaneous mea- 
surement of x and p on a single system gives a fidelity 
1/2, a well-known result. 

It is worth noting that optimally cloning squeezed 
states requires a variant of these SGCs, just as in Jl2[| . Let 
us consider for instance a family of quadrature squeezed 
states with squeezing parameter r. For such a family, the 
best symmetric doner must have the form of Eq. ([!]), but 
using the definition (3 = (^ + iap)/^/2 with a — exp(r). 



2 



These doners naturally generalize the SGCs and gives the 
same cloning fidelity, Eq. (||), for those squeezed states. 

In conclusion, we have established a link between op- 
timality of N — ► M symmetric Gaussian doners and the 
impossibility of simultaneously measuring two conjugate 
observables x and p. This results in a lower bound on the 
noise induced by cloning. The optimal cloning fidelity for 
coherent states was then derived, and was found to be in- 
dependent of which coherent state is to be copied. The 
optimal cloning of squeezed states was also found to be 
equivalent to that of coherent states, as expected since 
the former can always be obtained by applying a canon- 
ical transformation on the latter. It is unknown whether 
a doner specifically devised for other classes of states 
might yield a fidelity exceeding Eq. (|^). However, since 
minimum-uncertainty states are the closest to classical 
states, we conjecture that SGCs achieve the best possible 
fidelity if we require the doner to be covariant under ro- 
tations and translations in the phase space. Finally, even 
though the explicit transformation achieving the 1 — ► 2 
optimal SGC is known finding the N — > M cloning 
transformation that attains the maximum fidelity is still 
an open question. 

Appendix. We now prove that the variances of two cas- 
caded doners add. Consider a iV — ► M SGC, followed 
by a M — > L SGC. Let p be an arbitrary density op- 
erator supported on Tt® M . Since it is self-adjoint and 
compact, p has a denumerable spectrum: it can be ex- 
panded as p = <M£i)(&| with = <%, A 4 > 
and Ai = 1. Note that Ve > 0, 3d such that 
I X)i=i — 1| < e - Therefore, the output of the first 
doner can be decomposed as pu = Pd+ e d Bd where pd — 
J2i=i ■M&X&I i s supported on a d-dimensional subspace 
of Tt® M , Bd is a bounded operator, and lim^oo £d = 0. 
Since pu belongs to the symmetric subspace of TL® M , 
so will pd- Hence, we know that we can write pd in 
the form of a pseudo-mixture of pure product states 
Pd = J2i=i a i\4>t M )( ( t ) t M \ where the coefficients a« are 



not necessarily positive but satisfy ^2 i= 
or 0). Thus, when cloning a state 
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we have 
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Then, since the cloning map Cn,m is linear, cascad- 
ing the two doners yields Cm,lC n ,m{\^® N ){^® N \) = 
EiOiCM^f^i^D + edCMMBd)- As this expres- 
sion is a density operator (thus bounded) and the first 
term of its r.h.s. is positive, Cml(-B<j) must be bounded. 



Thus, the second term of the r. h. s. of Eq. ( |TT| ) becomes 
negligible when d — > oo. Now, using Eq.(|lj), we have 

Ttl-iCmxC N m{\^ N ){^ N \) = 
1 

n "M,L a N,M 



d 2 1 d 2 (3 e-lTl V^i-M /°3f.M 
xDfr + flMMDtfr + fl+Ofa) (12) 



with limd^oo rjd — 0. A little algebra then shows that 
this last expression is a Gaussian mixture centered on 
the original state whose variance is a 
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